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Abstract 



We observe hadronic W decays in ti — ► W(— ► + > 4 jet events using a 
109 pb _1 data sample of pp collisions at y/s = 1.8 TeV collected with the Collider 
Detector at Fermilab (CDF). A peak in the dijet invariant mass distribution is 
obtained that is consistent with W decay and inconsistent with the background 
prediction by 3.3<r. From this peak we measure the W mass to be 77.2 ± 4.6 
(stat+syst) GeV/c 2 . This result demonstrates the presence of two W bosons in ti 
candidates in the W{-^ > 4 jet channel. 



PACS numbers: 14.65.Ha, 13.38.Be 
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CDF presented the first direct evidence for top quark production with 19.3 pb -1 of data 
collected in 1992-93 |l[| and also reported a detailed kinematic study of top decays with 
this data sample ||. In 1994-95, the existence of the top quark was firmly established by 
CDF || and D0 [[|. CDF observed and studied it production using a sample of 67 pb -1 
employing techniques similar to those previously published [|]. [| . 

In the Standard Model, a top quark decays predominantly into a bottom quark and 
a W boson. As a part of the CDF kinematic studies of ti production and decay, we 
have searched for hadronic W decays in a ti -enriched sample containing a high-Pr lepton 
and four or more jets, which is consistent with the process pp —>■ tiX — > W + bW~bX — > 
lubjjbX with I = e or fi, and j representing a jet. Hereafter we call these events U W+ > 4 
jet events" . The observation of hadronic W decays in W+ > 4 jet events demonstrates 
the presence of two W bosons in the final state for ti candidates. 

The studies reported here use a 109 ± 7 pb -1 data sample. Three analysis techniques 
yield signal-to-background ratios ranging from 0.23 to 1.4 in the dijet invariant mass 
region 60-100 GeV/c 2 . 

The CDF detector consists of a magnetic spectrometer surrounded by calorimeters 
and muon chambers @. A silicon vertex detector (SVX), located immediately outside 
the beampipe, provides precise track reconstruction in the plane transverse to the beam 
and is used to identify secondary vertices from b and c quark decays ||. The momenta of 
charged particles are measured in the central tracking chamber (CTC) which is in a 1.4 T 
superconducting solenoidal magnet. Outside the CTC, electromagnetic and hadronic 
calorimeters cover the pseudorapidity || region \t)\ < 4.2 and are used to identify jets 
and electron candidates. Outside the calorimeters, drift chambers in the region \r)\ < 1.0 
provide muon identification. 

To select ti — > W+ > 4 jet candidates, we require an isolated, high-P T electron or 
muon (P T > 20 GeV/c), high missing transverse energy || (]fi T > 20 GeV), three jets 
with E T > 15 GeV and \rj\ < 2.0, and a fourth jet with E T > 8 GeV and \rf\ < 2.4. After 
these cuts 163 events remain. The expected fraction of ti events in this sample is 33%. 
The fraction was determined by extrapolating the background calculation performed on 
the W+ > 3 jet sample for the ti cross section measurement to include the additional 
requirement of a fourth jet. 

For the event selection, neither the jet Ex nor the Ifix are corrected for detector ef- 
fects. However, when calculating a dijet invariant mass, jet energies are corrected by a 
pseudorapidity and energy- dependent factor which accounts for such effects as calorimeter 
nonlinearity, reduced response at detector boundaries, contributions from the underlying 



event and multiple interactions and losses outside the clustering cone [11]. We apply an 



additional energy correction to the four highest-P^ jets in a it candidate event. The cor- 
rection depends on the type of parton they are assigned to: a light quark, a hadronically 
decaying b quark, or a b quark that decayed semileptonically |l|. This parton-specific cor- 
rection was derived from a study of ti events generated with the HERWIG Monte Carlo 
program []TJ] . This correction is applied in order to reconstruct the original parton energy 
as closely as possible. These cuts and corrections are the same as those used in the top 
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mass measurement U before requiring a b jet. 

Monte Carlo it events are generated by the HERWIG program ]12| with a top mass 
of 175 GeV/c 2 |0. The expected background from direct W+ jets production is esti- 



mated using the VECBOS Monte Carlo calculation [14| to generate W+ 3 parton matrix 
elements. To model W+ > 4 jet events we implement a simulation of parton fragmen- 
tation with a shower algorithm based on HERWIG. The lowest-order matrix elements 
are sensitive to the choice of the mass scale in the strong coupling constant a s . We use 
two Q 2 scales, namely, the square of the average jet Pt ((Pt) 2 ) and the square of the W 
boson mass (M^). We use {Pt) 2 as a standard Q 2 scale and for an estimate of the 
systematic uncertainty. A more detailed description of the Monte Carlo samples can be 
found elsewhere IJ. 

The lepton-^r transverse mass distribution is shown in Figure |l[ Also shown are the 
Monte Carlo expectations from it as well as it plus QCD W + jets background. The 
Monte Carlo expectations are normalized to the number of events observed. 

We use three methods to search for W decay to two jets: a top mass reconstruction 
technique, a total transverse energy cut, and the identification of both b jets. 

In the first method, we search for a W decaying to two jets by kinematically recon- 
structing it — > ivbjjb with a constrained fitting technique similar to the one used for the 
determination of the top mass PJ. The constraints used in the top mass measurement 
include the requirement that the two jets hypothesized to come from the W decay have 
an invariant mass equal to the W mass. Here we eliminate the W mass constraint. There 
are multiple solutions due to the ambiguity in determining the neutrino longitudinal mo- 
mentum and the assignment of jets to the parent partons. We choose the solution with 
the lowest \ 2 . The resolution and signal-to-background ratio for tt events is improved by 
requiring a 6-tag in the event and requiring that the tagged jet correspond to a b jet in the 
fit. To tag b jets, we exploit either the long lifetime of b quarks by requiring a secondary 
vertex (SVX 6-tag), or the semileptonic decays of b quarks by searching for additional 
leptons (SLT 6-tag) The fraction of correct W jet assignments was found from Monte 
Carlo studies to be improved from 22% before 6-tagging to 37% after tagging. 

Of the 163 events passing the W+ > 4 jet selection cuts, 37 have at least one 6-tag found 
by the SVX or SLT algorithms. The expected background is 8.9^;° events, calculated 
using techniques described in JIIJ. The dijet invariant mass distribution for these events is 
shown in Figure Also shown are the Monte Carlo distributions for tt and background, 
normalized to 28.1 and 8.9 events respectively. We note a strong enhancement in the data 
distribution near the world-average W mass Jl5 . 



The second method employs a cut on the minimum total transverse energy Ht to 
enhance the tt contribution in the sample. Ht is defined as the scalar sum of the Pt of 
the lepton (E T for an electron), the four highest-^ jets, and the IpT @- We optimize 
the Ht threshold (310 GeV) using tt and background simulations. This cut improves the 
fraction of tt events in the sample from 33% to 55%. 

We calculate the invariant mass of each two-jet pair out of the four highest--Er jets. 
Monte Carlo simulations predict that two jets from W decay are included in the four 
highest- .Et jets with an efficiency of 66%. There are six combinations of two jets in 



7 



each event. To determine the contributions from non-tt backgrounds and combinatoric 
backgrounds in ti events, we fit the dijet mass distribution Mjj outside the W mass 
region to the following function: N bg f bg (Mjj) + N comh f comb (M j:j ), where N hg and N comb 
are the number of non-tt and ti combinatoric background dijet pairs respectively, and 
f bg and fcomb are their dijet mass distributions, derived from Monte Carlo. Since non-ti 
backgrounds, such as WW, WZ, and bb are smaller and give dijet mass distributions 
similar to QCD W + jets, f bg is obtained from VECBOS calculation of QCD W + jets 
background. 

We extract the W — > 2 jets signal by subtracting the QCD W+ jets background and 
the combinatoric background as shown in Figure ^. For the W+ > 4 jet sample of ti 
production, the reconstructed mass peak of the W decaying into dijets is well described by 
a Gaussian distribution with a standard deviation of 11.7 GeV/c 2 . We fit this distribution 
to a Gaussian function with a fixed width and obtain 29 ±13 W — > 2 jet events in the mass 
region 70-90 GeV/c 2 , which is consistent with the number of W — > 2 jet events (20 ± 6) 
expected from the CDF ti production cross section measurement [l(|. This excess is 



inconsistent with the expected background by 2.8a, corresponding to a probability of 2.6 
x 10~ 3 of it being a background fluctuation. The signal-to-background ratio in the mass 
region 60-100 GeV/c 2 is 0.23, and the overall fraction of non-ti background is (27 ± 14)%. 
The dijet invariant mass peak has a mean of 77.1 ± 3.8(stat) ± 3.6(syst) GeV/c 2 . The 
systematic uncertainty comes mainly from the jet energy scale, which is the uncertainty 
on how well our measured jet energies correspond to the original parton energies. The 
systematics are listed in Table [I], and are described in more detail in Ref. |L3 . 



In the third method, we isolate the hadronic W decay by identifying two b jets in 
the four highest-E'T jets of W+ > 4 jet events. Such events are hereafter called "double 
6-tagged" W+ > 4 jet events. Double 6-tagging eliminates ambiguity about which two 
jets are assigned to the W, and further suppresses non-tt background. Once we find at 
least one b jet with an SVX tag or an SLT tag in an event, we look for the second b jet. 
The second 6-tag can be an SVX tag, an SLT tag or a Jet Probability tag The Jet 
Probability algorithm obtains the probability that a jet is consistent with the decay of a 
zero-lifetime particle by using the impact parameters of the tracks in the jet as measured 
in the silicon vertex detector. Jets with probabilities less than 5% are considered tagged. 
We form the invariant mass of the remaining two untagged jets. 

In our HERWIG ti Monte Carlo sample, 25% of the events are double 6-tagged. The 
technique finds the two correct W jets in 43% of these events. However, this number 
is sensitive to the amount of initial and final state radiation in the simulation, since 
the largest contamination (38% of double 6-tagged events) comes from events where the 
four highest-.&r jets do not correspond to two b jets and two jets from W decay. Other 
ti backgrounds include events where both W bosons decayed leptonically and only one 
lepton was identified (8%), events where a c quark from W decay was tagged (8%), and 
mistags in jets from W decay (3%), where a mistag means a tag of a jet originating from 
neither a b nor a c quark. 

We expect 1.3 ± 0.3 double-tagged events from non-tt background, with the largest 
contributions being 0.6 ± 0.2 events from mistags and 0.4 ± 0.2 events from Wbb and 
Wcc processes which are QCD W+ jets events containing real heavy flavor. 
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In the data, we find eleven double 6-tagged events. The dijet invariant mass of the 
two untagged jets in these eleven events is shown in Figure Eight of the eleven dijet 
combinations fall in the mass window of 60-100 GeV/c 2 . In Monte Carlo simulations of 
both QCD W + jets and it backgrounds, only about a third of the dijet mass combinations 
fall in this window. The inset plot in Figure |] shows the mass of the hadronic W candidates 
against the transverse mass of the leptonic W candidates in the same events. 

We fit this mass distribution to a sum of a Gaussian W — > 2 jets signal, ti backgrounds, 
and non-ii backgrounds, and obtain a W mass of 78.1 ± 4.4(stat) ± 2.9(syst) GeV/c 2 , 
with the systematic uncertainties listed in Table [I]. The fitted W signal has 8.7 events, 
and 0.7 events from ti backgrounds. Constraining the fit to the expected W — > 2 jets 
fraction in ti events yields 6.4 W events, and a W mass of 78.3 ± 5.1(stat) ± 2.9(syst) 
GeV/c 2 . We use this value for obtaining a final combined W mass. 

We use likelihood fits to evaluate the significance of this peak. We first fit the data 
using a Gaussian W signal whose mean is the world-average W mass plus ti and non-it 
backgrounds. Next we remove the signal term and note the change in likelihood. We 
then repeat this procedure on Monte Carlo pseudo-experiments using only ti and non-it 
background events, and study how often the likelihood changes by at least the amount 
seen in the data. This change in likelihood was seen with a probability of 1.7 x 10~ 3 , 
corresponding to a Gaussian significance of 2.9a. If we construct the pseudo-experiments 
with the expected W signal fraction, we find that a change in likelihood larger than the 
one seen in the data occurs 15% of the time. 



Uncertainty (GeV/c 2 ) cut double 6-tag common combined 

Statistical 4.6 5.1 - 3.5 

Jet energy scale 2.8 

(a) Detector effects 1.9 1.9 1.9 

(b) Soft gluon effects 2.0 2.0 2.0 

Other systematics 0.6 

(a) Backgrounds 0.5 0.2 

(b) Fitting 1.9 0.6 

(c) H T cut 1.0 

Total uncertainty 5.8 5.9 2.8 4.6 



Table 1: The list of uncertainties (in GeV/c 2 ) in the W — > 2 jets mass for the 
cut and double 6-tag analyses. The columns "common" and "combined" correspond to 
the uncertainties common in the two analyses and the combined uncertainties of the two 
analyses, respectively. 

An overall result is obtained by excluding the eleven double 6-tagged events from the 
Ht cut analysis. We then combine the two analyses to obtain a W mass peak significance 
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of 3.3cr, corresponding to a probability of 5.4 x 10 -4 of it being a background fluctuation 
and a W mass of 77.2 ± 3.5(stat) ± 2.9(syst) GeV/c 2 with the systematic uncertainties 
listed in Table |l|. This is consistent with the current world-average W mass of 80.375 ± 
0.120 GeV/c 2 0. 

In conclusion, we observe hadronic W decays as a dijet mass peak in ti — > W+ > 4 
jet events. This demonstrates the presence of two W bosons in the ti candidates. In the 
future, these techniques can be used to set limits on nonstandard top decays. 

We thank the Fermilab staff and the technical staffs of the participating institutions 
for their contributions. This work was supported by the U.S. Department of Energy 
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the Ministry of Science, Culture, and Education of Japan; the Natural Sciences and 
Engineering Research Council of Canada; the National Science Council of the Republic 
of China and the A. P. Sloan Foundation. 
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Figure 1: The lepton-^y (neutrino) transverse mass distribution in the W+ > 4 jet data 
(points) along with the Monte Carlo expectations from ti (dotted) and tt plus QCD W 
+ jets background (dashed). 
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Figure 2: Invariant mass distribution for the two jets that are associated with the W 
by the constrained kinematic fit. The shaded histogram shows the expected distribution 
from background events; tt plus background events are shown as a dashed line. 
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Figure 3: Dijet mass distribution in W+ > 4 jet events after the Ht > 310 GeV cut. 
Also shown are a fitted curve of the sum (solid) of the VECBOS W + jets background 
(shaded) and the combinatoric background from ti events. A Gaussian distribution fitted 
to the W mass peak is shown by the dotted curve. The inset shows the W mass peak 
after subtracting the backgrounds. 
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Figure 4: Dijet mass spectrum of the two untagged jets in double 6-tagged events. The 
shaded curve shows the expected distribution from QCD W + jets background and the 
dashed curve is from tt plus background. The inset compares the lepton-^T (neutrino) 
transverse mass to the dijet invariant mass in these eleven events showing evidence for 
events with both a leptonically-decaying W and a hadronically-decaying W. 
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